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ABSTRACT 

Vorticity forms the basis for a new approach to the anal- 
ysis of inlet distortion in axial flow fans and compressors. 
Farmer first formulated this approach and this paper repre- 
sents a test of the usefulness of his theory. A summary of 
recent developments in the field along with a discussion of 
the effects of distortion on stall margin forms a background 
for this complex problem. A computer program calculates the 
vorticity pattern at the compressor face using data that are 
read from a magnetic tape. The data, which were provided by 
NASA Lewis Research Center in Cleveland, consist of eleven 
stall events. The output of the program contains tables of 
stagnation pressure and the partial derivatives of pressure 
in the R and theta directions as well as three maps. The 
three maps are of pressure and radial and circumferential 
vorticity. The results of the program show a correlation be- 
tween a ring of large positive circumferential vorticity and 
stall, leading to the conclusion that the stall was caused by 
the increase in blade loading induced by the vorticity. This 
conclusion suggests a formulation for a universal inlet 
distortion index and provides a basis for the evaluation of 
the vorticity approach to the problem. 



2 



TABLE OF CONTENTS 



I. INTRODUCTION — 8 

II. NATURE OF PROBLEM 10 

A. SOME PERTINENT DEFINITIONS * 10 

B. RECENT DEVELOPMENTS IN THE FIELD 14 

C. VORTICITY APPROACH 20 

III. INTERPRETATION OF VORTICITY MAPS 25 

A. EMPIRICAL APPROACH 25 

B. ANALYTICAL APPROACH 33 

IV. A DISTORTION INDEX 37 

A. DEVELOPMENT OF THE INDEX 37 

B. TESTS OF THE DISTORTION FACTOR 40 

V. SUMMARY AND REVIEW 46 

VI. CONCLUSIONS 47 

APPENDIX A EXPLANATION OF DATA 49 

APPENDIX B ACCURACY OF COMPUTER PROGRAM 55 

COMPUTER PROGRAM — : 65 

LIST OF REFERENCES 84 

INITIAL DISTRIBUTION LIST 86 

FORM DD 1473 93 



3 



LIST OF TABLES 



I Distortion Factors 19 

II Possible Combinations in a Rotor 34 

A- 1 Steady-state Operating Conditions 53 

A-II Times of Interest 54 

B-I Predicted Values of a Function vs. Actual Values-- 57 

B-II Test of Print Routine Symbols 53 



4 



LIST OF FIGURES 



1. Compressor Map, J-85-GE-13 

2. Compressor Stability 

3. Instantaneous Pressure Map 

4. Ambiguous Direction of Axial Vorticity 

5. Computer Program Flow Diagram 

6. Instantaneous Radial Vorticity Map 

7. Scale for Vorticity Maps 

8. Case I Circumferential Vorticity 

9. Case II Circumferential Vorticity 

10. Case III Circumferential Vorticity 

11. Case IV Circumferential Vorticity — 

12. Equivalence of Blade Circulation and Vorticity 

13. Additional Lift Due to Radial Vorticity 

14. Vorticity Pattern 

15. Induced Downwash 

16. Regions of Compressor Face 

17. Distortion Index, Case I 

18. Distortion Index, Case II 

19. Distortion Index, Case III 

20. Distortion Index, Case IV 

21. Vorticity Pattern Produces Large Index, No Stall 

A-l Schematic of Digitizing System 

A-2 Placement of Probes 

B-l Pressure Vs. Theta, r'=0.41 

B-2 Pressure Vs. Theta, r'=0.59 



11 

12 

17 

21 

23 

27 

28 

29 

30 

31 

32 

34 

34 

35 

36 

36 

42 

43 

44 

45 

41 

50 

51 

58 

59 



B-3 Pressure vs. Theta, r'=0.73 60 

B-4 Pressure vs. Theta, r'=0.85 61 

B-5 Pressure vs. Theta, r'mO.95 62 

B-6 Map of Circumferential Vorticity 64 



6 



ACKNOWLEDGEMENTS 



The author wishes to express sincere appreciation for 
the patience and guidance of Professor Allen E. Fuhs of the 
Department of Aeronautics at the Naval Postgraduate School 
during the completion of this thesis. The author also is 
grateful for the assistance of Hans Doleman of the W. R. 

Church Computer Center at the Naval Postgraduate School. 

Several persons and organizations provided data and/or 
assistance for this project. The author wishes to thank 
Captain Barry Browns tein of AFAPL/WPAFB for his contributions 
to the development of the computer program. The' author also 
expresses appreciation for the assistance of James Calogeras, 
Robert Coltrin, Paul Burdstadt and members of the staff at 
NASA Lewis Research Center in Cleveland, Ohio and for their 
willingness to supply the data which were used in this thesis. 



7 



I. 



INTRODUCTION 



The problem of surge in axial flow fans and compressors 
in jet engines is a serious one. Unsolved it can result in 
the loss of a multi-million dollar aircraft and possibly in 
the loss of the life of the pilot. Manufacturers are there- 
fore forced to design compressors which operate with decreased 
performance. 

During the early years of the use of axial compressors in 
aircraft jet engines, before high speed flight was a common 
occurrence, the problem of surge was treated as a steady-state 
problem with satisfactory results. However, with the advent 
of high performance and high speed military aircraft, and in 
particular the F-lll which employed the TF-30 turbo- fan engine, 
[Ref. 1] it became evident that the problem could no longer be 
considered to be steady-state. There were some attempts made 
to . improve the steady-state distortion factor, but these were 
unsuccessful. During the next few years, various attempts 
were made to develop some sort of a distortion index using the 
fluctuations of the stagnation pressure as a starting point. 
These attempts met with some degree of success, but none was 
the complete answer to the problem. 

It has been proposed that a distortion factor which uses 
vorticity as a measure of instantaneous distortion would 
provide insight into the fluid mechanics of the problem. 

Farmer [Ref. 2] first suggested this approach and developed 
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approximate equations for the vorticity at the compressor 
face in terms of measured quantities, i.e., stagnation 
pressure data. 

Since most of the distortion factors developed by other 
researchers were the outgrowth of instantaneous maps of the 
pressure distribution over the compressor face, it seems 
logical to take that same route to the development of a dis- 
tortion index based on vorticity. Therefore a computer pro- 
gram was developed to calculate the instantaneous vorticity 
pattern and produce the maps of the pattern on the printer 
component of the computer. This paper will attempt to inter- 
pret the output of this program and prove the validity of 
this approach. Suggestions then will be made concerning the 
formation of a distortion factor that will allow the pre- 
diction of surge for a wide variety of engines. 
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II. NATURE OF PROBLEM 



A. SOME PERTINENT DEFINITIONS 

In order to understand the complexity of the problem, one 
must first be familiar with some of the nomenclature and basic 
concepts that are used in the field of compressor stability. 
Brimelow [Ref. 3] gives an excellent presentation of the 
present theories on this subject as well as a concise 
explanation of the nomenclature. 

The most likely place for an engine to become unstable is 
in the fan or compressor. These two components are said to 
be in surge when several or all of the blades become too 
heavily loaded and are forced to operate in an aerodynamically 
stalled condition. If one examines the compressor map of the 
J-35-GE-13 jet engine in Fig. 1, he can see that the conditions 
for surge are represented by a line on the compressor map. 

This line is generally called the stall or surge line. The 
most efficient operating condition of the compressor is as 
close to this stall line as possible without stalling. In 
order to insure that the compressor will not surge, it must be 
run with a safety margin between the operating point and the 
stall line. The vertical difference between the operating 
point and the stall line is defined as the stall or surge 
margin. In today's high performance aircraft it is necessary 
to have as small a surge margin as possible in order to de- 
crease the specific fuel consumption and increase the total 
thrust, while building the engine to be as light as possible. 
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I 

FIG. 1 COMPRESSOR MAP J-85-GE-13 



It appears to be an easy task to define a small surge mar- 
gin and then run the engine at a point such that this surge 
line will never be crossed. It turns out, however, that the 
compressor map is only a representation of the steady-state 
operating conditions of the compressor. Once that fact is 
known, it is understandable that fluctuations in the flow may 
affect the stall margin and cause the compressor* to stall even 
though the steady-state operating point may have been well 
below the stall line. 

The term compressor instability may best be defined through 
visual means. Figure 2 is a hypothetical compressor map with 
only one constant speed line. On this line six operating 
points are labeled Pi, P2, Rl , P3, P4, and R2 . If the com- 
pressor is operating on the negative slope of this line at 
point Rl, its operation is stable. It is slable because the 
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operating point will tend to return to R1 if it is forced in 
either direction along the constant speed line. If it is 
forced to P2 the mass flow will be decreased. The character- 
istics of the chamber behind the compressor, which is known 
as the receiver, require a lower pressure ratio at the re- 
duced mass flow rate. The compressor discharge flow will be 
forced to expand into the receiver. This will cause an in- 
crease in mass flow rate and the operating point will return 
to Rl. In a similar manner if the operating point is forced 
to PI, it will be forced to return to Rl. In this case the 
pressure ratio will be less than that required by the receiver 
at the increased mass flow rate and a back pressure will re- 
sult. The back pressure will cause a reduction in flow rate 
and the operating point will return to Rl. 




FIG. 2 COMPRESSOR STABILITY 
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If the compressor is operating on the positive slope of 
the constant speed line at point R2 , the operation of the 
compressor is unstable. In this case, if the operating point 
is disturbed from R2 , it will continue to move in the direc- 
tion in which it initially was moving and will not return to 
R2 . If the operating conditions move to point P4, a decrease 
in pressure ratio will occur. This will in turn set up a 
back pressure since the pressure in the receiver will be 
greater than the compressor discharge pressure, and the flow 
rate will be reduced. The decrease in flow rate will con- 
tinue to reduce the pressure ratio, and the operating point 
will continue to move away from R2 . If the operating point 
moves in the other direction along the constant speed line, 
the pressure ratio will increase which, in turn, will in- 
crease the flow rate, and again the operating point will 
continue to move away from the design point R2 . 

One can see, therefore, that the limit for stable com- 
pressor operation is at the point where the slope of the 
constant speed line on the compressor map goes to zero. The 
compressor must operate where the slope of the constant speed 
line is negative and not be allowed to reach the pressure 
ratio that corresponds to the limiting point for stable com- 
pressor operation. 

With high performance aircraft the transient distortions 
induced by sources external to the engine and in the engine 
inlet will be of a large magnitude and, therefore, more 
dangerous to the stability of the compressor. The turbulence 
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is likely to be generated by some or all of several sources. 
Some of these sources are: shock waves, boundary layers, hot 

gases emitted from rockets launched from the aircraft, and 
off-design operation of the inlet. Since the stall margin is 
such a critical variable in the operation of a jet engine, it 
is not feasible merely to increase the stall margin to such 
an extent that time-variant distortion cannot cause the 
compressor to stall without a disastrous drop in the overall 
efficiency of the engine. It is necessary to obtain more 
information on the source of the distortion, what effect the 
distortion has on the inlet-compressor interface, and just 
exactly what type of distortion is most likely to cause a 
compressor to surge. 

B. RECENT DEVELOPMENTS IN THE FIELD 

With the advent of the new high performance aircraft it 
became evident that the method that was used for the evaluation 
of the inlet-compressor interface was no longer effective. At 
that time this evaluation had been done wholly by an analysis 
of the steady-state distortion through the use of a steady- 
state distortion factor. In a high performance aircraft the 
inlet flow is more likely to be turbulent making it necessary 
to evaluate the effects of time variant distortion on the 
compressor. Since a turbulent flow is a function of four 
variables, which in a cylindrical coordinate system are R, 9 , 

Z and time, the problem became a much more complex one. The 
magnitude of the fluctuations in pressure and velocity in a 
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turbulent flow is such that they cannot be considered small 
with respect to the steady- state values. 

In the light of turbulent inlet flows, it becomes apparent 
why the steady-state distortion index was not adequate. It 
does not take the time variable into account. When the inlet 
flow of a compressor is turbulent, the -steady-state distortion 
factor can be well within the tolerences for stall, and yet 
the compressor may stall due to the random velocity and 
pressure fluctuations in the flow. 

The first attempt to solve this problem resulted in an im- 
proved steady-state distortion factor. As might be expected, 
however, that factor was also inadequate for the solution of 
the problem. The amount of turbulence is often quantified 
by the variation of root-mean-square stagnation pressure 
(A ' Still reluctant to tackle the turbulent flow 

problem, researchers attempted to devise a method by which 
the improved steady-state distortion factor could be used 
with a correlation for the amount of turbulence [Ref. 1]. 

This method was also ineffective, but it gave some insight 
into the complexity of the problem. 

At this point, it was evident that some method of predict- 
ing the effects of time variant distortion must be developed 
to be used in the inlet-compressor matching problem. The 
method used was to introduce a system of instrumentation that 
had a quick response and was able to measure the stagnation 
pressure at any particular point in time. Measurements then 
were taken at small intervals of time with pirobes located on 



15 



pressure rakes at the compressor face. The pressures at 
these discrete points were converted to a distribution and 
presented in the form of an instantaneous pressure map. By 
fixing the variable Z at the compressor face, researchers 
eliminated one variable; and by taking instantaneous pressure 
readings, they made the problem two-dimensional at any point 
in time. The influence of time was expressed in the analysis 
by the time continuity of a set of maps. It became possible 
to find an empirical correlation between the instantaneous 
distortion, as represented by the instantaneous pressure map, 
and surge. In Ref. 5 Calogeras, Burdstadt, and Col trill found 
an empirical correlation between the instantaneous pressure 
distribution and the onset of surge. One of the maps they 
used in the development of the correlation is depicted in 
Fig. 3. Using their correlation between the instantaneous 
pressure distribution and surge, they developed an instan- 
taneous inlet distortion factor with which it was possible to 
predict stall with some measure of success. This factor was 
developed for the J-85-GE-13 engine with a particular inlet 
geometry and was not tested with either other engines or 
other inlets. It was tested at angle of attack and proved 
quite successful in that context. 

Plourde and Brimlowe [Ref. 6] also developed an instan- 
taneous inlet distortion factor. This factor was used for 
the TF-30 turbofan engine which is installed in the F-lll. It 
is interesting to note that both of these factors, for the J-85 
and the TF-30, are primarily circumferential distortion factors. 
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Another interesting point is that while each of these factors 
is useful in the analysis of the case for which it was devel- 
oped, neither would be useful, except in a general way, to 
anyone working on the inlet distortion problem in a different 
engine . 

The period of time for which a particular distortion 
pattern must persist in order to cause stall was found to 
be extremely small. Reference 7 states that for the TF-30 
this time may be approximately three milliseconds, while 
Ref. 5 has stated that the time scale for the J-85 can be as 
short as one millisecond. Reference 8 contains a method of 
calculating the length of time for which a distortion must 
persist in order to cause steady-state stall. In general 
this time, depending on the size of the engine, engine speed, 
and number of blades, is between one and five milliseconds. 

If the duration of a distortion is less than this time, stall 
will not occur, due to dynamic suppression of stall, even 
though the type of distortion might have caused a stall had 
it persisted for a longer period of time. Unfortunately the 
exact time that a distortion must persist cannot be calculated 
since it may be a function of many more variables than those 
mentioned above. 

Upon examination of several other papers, [9, 10, 11, 12] 
it becomes evident that the large amount of work that has been 
done on this problem has resulted in two conclusions. One of 
these results is an empirical relationship between instanta- 
neous inlet distortion and stall, while the other result is 
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lack of agreement on the quantitative form of the relation- 
ship. As a consequence dozens of instantaneous distortion 
factors are used at the present time. This lack of consis- 
tency can be realized by examining Table I which contains 
six distortion factors f all of which have been employed for 
the same engine at various points in time. 

TABLE I - DISTORTION FACTORS 
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Brownstein [Ref. 13] puts forth a possible universal in- 
let distortion factor using pressure fluctuations as a base 
for the index. It appears to have solved some of the prob- 
lems encountered by earlier factors and is a more general in- 
dex. This index seems to be very promising. Another, 
different approach to the problem was put forth by Farmer 
[Ref. 2] and again in Ref. 8. 

C. VORTICITY APPROACH 

The approach presented in Ref. 2 and Ref. 8 is new. 

Using Crocco 1 s theorem it was possible to obtain an ex- 
pression for vorticity at the compressor face. Then, with 
an order of magnitude analysis, using data taken for the 
J-85-GE-13 jet engine. Farmer was able to eliminate some of 
the terms involved in the original equation and emerge with 
a simple, easily calculated expression for vorticity at the 
compressor face. The two expressions, for radial and 
circumferential vorticity are 



For a more detailed explanation, consult Farmer's work in 
Ref. 1. Due’ to the fact that vorticity lines do not end, 
i . e . , V X UJ = 0, the axial vorticity can be calculated vising 
a control volume. Unfortunately the direction of the axial 




( 1 ) 




( 2 ) 
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vorticity will be ambiguous. As shown in Fig. 4, it is not 
possible to determine if the axial vorticity enters the 
control volume from the front or the back. 




FIG, 4 AMBIGUOUS DIRECTION OF AXIAL VORTICITY 



Reference 2 states that circumferential vorticity is 
likely to have a larger effect on the stall margin than is 
radial vorticity. This is in concurrence with other work 
done on this problem, but it may be beneficial to consider 
the effects of radial vorticity in the light of the fact 
that the duration of a distortion pattern necessary to cause 
stall is extremely short. It should be noted, however, that 
radial vorticity in a rotor is an unsteady phenomenon [Ref. 8] 
and therefore very difficult to treat. 

This new approach seems to have merit, and its usefulness 
should be tested using some of the large supply of existing 
data that have been accumulated in the study of this problem. 
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It would appear that some method of producing instantaneous 
vorticity maps similar to instantaneous pressure maps would 
be useful in the evaluation of this theory. Also should the 
theory be confirmed by these maps, the maps would be instru- 
mental in the formation of a distortion index based on 
vorticity. 

The high speed digital computer is the best choice for 
the means of production of these maps. The main .reason for 
this choice is the amount of time that would be necessary 
for their production by hand or even on the most sophisti- 
cated desk calculator that is available today. Also, it is 
necessary that the partial derivatives of pressure in the 
r’ and € directions be calculated, and this can be done most 
accurately by the use of numerical curve fitting routines 
that already have been written for the computer. 

The problem then is to develop a computer program that 
will fulfill the following requirements: read- the data from 

a magnetic tape, define a function in terms of r’ and 0 that 
will provide the values of the normalized pressure and its 
two derivatives at any point on the compressor face, calcu- 
late the normalized axial velocity, calculate the radial and 
circumferential vorticity at any point on the compressor 
face, and finally print maps of pressure and radial and 
circumferential vorticity at the compressor face. These re- 
quirements are shown more clearly by the flow chart in Fig. 5 

Once the program lias been developed and meets the require 
ments mentioned above, the maps should be interpreted as to 
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the effect of the instantaneous vorticity distribution on the 
stall margin. It is hoped that there will be some relation- 
ship between the instantaneous vorticity pattern and surge 
other than the purely empirical correlation between the 
instantaneous pressure distribution and surge that has been 
found before. 




FIG, 5 COMPUTER PROGRAM FLOW DIAGRAM 

Farmer, Iverson and Fuhs [Refs. 8 and 14] have presented 
one such logical relationship between vorticity and stall. 

It seems that vorticity can cause increased blade loading. 
Obviously, if the increased blade loading is of a large 
magnitude and has a sufficiently long duration, this vorticity 
could be the cause of compressor stall. 

If it is not possible to find a logical correlation be- 
tween instantaneous vorticity and stall, it still might be 
possible to find an empirical correlation between vorticity 
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and stall more satisfying than the relationship that has been 
discovered between instantaneous pressure and stall. This 
would seem to be a logical expectation since vorticity con- 
tains terms of a higher order than pressure and may be a 
truer measure of distortion. 

These two approaches to the interpretation of the 
vorticity maps both have merit, and both should be considered 
in this paper. 

The discussion of this topic will contain three more 
considerations. First, it should contain some type of 
quantitative evaluation of the utility of the instantaneous 
vorticity distribution in the analysis of inlet distortion. 
Secondly if the vorticity map proves to be a useful tool, 
suggestions as to the most beneficial employment of the map 
in the further study of the problem will be made, and some 
work will be done on the development of an inlet distortion 
factor. Finally some insight will be provided into possible 
alternative paths to the solution of the problem. 
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III. INTERPRETATION OF VORTICITY MAPS 



A. EMPIRICAL APPROACH 

A computer program which produces maps of instantaneous 
vorticity at the compressor face lends itself quite readily 
to the empirical approach to the interpretation of these 
maps. It produces a large number of maps that are easily 
examined. If a pattern of vorticity is present in several 
of the events just prior to stall, and not present at other 
times, it is likely that the pattern represents the distortion 
which was the cause of compressor surge. In the examination 
of these maps it is necessary to consider the duration of the 
pattern which seems to be causing surge. If the pattern does 
not persist for a sufficient length of time, its effects on 
the compressor are likely to be minimized by dynamic suppres- 
sion of stall. It is possible to compute a number, known as 
the reduced frequency, which will provide some insight into 
the influence of dynamic suppression on a distortion pattern's 
ability to cause compressor surge. Equation (3) , from Ref. 

14, 

Reduced Frequency = .J 1-JL9 c = chord (3) 

(J = tip speed 

is the equation for reduced frequency. For the J-85 the 
length of time that a pattern must exist in order to be 
relatively free from the influence of dynamic suppression of 
stall is approximately one millisecond. This number, however, 
is more of an order of magnitude approximation than a hard 
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and fast rule. It could be that a particularly strong 
distortion pattern with a duration of slightly less than 
one millisecond would be quite capable of inducing stall in 
the compressor. For a more complete explanation of dynamic 
suppression of stall, consult Carta [15], Erikson and Redding 
[16] and Carta [17], 

Examination of the maps generated by the computer program 
leads to the conclusion that the radial vorticity pattern is 
unlikely to yield an empirical correlation to stall or least 
not one that is evident through visual inspection. The 
pattern of radial vorticity is generally quite random, and 
the same pattern does not persist for more than a fraction 
of a millisecond. Figure 6 contains one of the radial 
vorticity maps drawn by the program. Figure 8 contains the 
scale to be used in examination of the maps. The numerical 
values in Fig. 7 are the values of the constant normalized 
vorticity lines that are represented by the boundary between 
the two symbols on either side of the number. This scale is 
the same for all maps of vorticity that are presented in 
this thesis. 

The patterns of instantaneous circumferential vorticity 
prove much more enlightening. The maps in Figs. 8-11 are 
patterns that appear prior to stall in the first four of the 
eleven stall events. 

The ring of large positive vorticity at the tip of the 
blades and extending nearly an entire 360 degrees seems to 
be quite likely to have caused stall. It is not feasible to 
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FIG, 9 CASE II CIRCUMFERENTIAL VORTICITY 
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examine every map in a stall event, due to the amount of 
computer time required for the production of these maps; 
each set of three maps requires 50 seconds of computer time. 
With 2048 sets of three maps per stall event, this would 
amount to 28 hours of computer time for each event. A 
random sampling of vorticity maps was taken, and the ring 
pattern does not appear except at such a time that it could 
be the cause of stall. The similarities between these maps 
are striking. A distortion index based on this pattern 
could yield some insight into how the compressor stalls. 

It would appear from this pattern that the tips of the blades 
are the cirtical point of the compressor for the J-85. 

B. ANALYTICAL APPROACH [Ref. 8] 

There has been a large amount of work done for predicting 
the effects of nonuniform flows on the performance of a 
cascade. Reference 8 contains an excellent summary of some 
of the results of this work as well as sample calculations 
of the effects of secondary flow on blade loading. Further 
insight into these effects can be gained by consulting Ref. 14. 

The primary effects of radial vorticity can be understood 
easily. Reference 18 states that an airfoil can be represened 
by a reentrant vortex. It is fairly easy to see that the 
blades in a compressor could be approximated by a group of 
negative radial vortices. (See Fig. 12). It is then a logical 
supposition that positive radial vorticity would tend to de- 
crease the blade loading and negative radial vorticity would 
tend to increase blade loading , as in Fig. 1.3. The trouble 
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with this observation is that radial vorticity in a rotor 
is an unsteady phenomenon. Table II classified types of 
vorticity in a rotor and stator as steady or unsteady. 




FIG. 12 EQUIVALENCE OF BLADE CIRCULATION AND VORTICITY 





FIG. 13 ADDITIONAL LIFT DUE TO RADIAL VORTICITY 



TABLE II. POSSIBLE COMBINATIONS IN A ROTOR [Ref. 8] 



Vortex 


Direction pf Vortex 


Flow 


Sheet 


Radial 


Unsteady 


Filament 


Radial 


Unsteady 


Filament 


Circumferential 


Unsteady 


Sheet 


Circumferential 


S teadv 
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The influence of circumferential vorticity on blade 
loading is not so clear. In Ref. 8 the secondary effects 
of one vorticity pattern were calculated with the following 
results. It was found that circumferential vorticity could 
cause additional blade loading through its secondary flow 
effects. The vorticity pattern for which these effects were 
calculated is remarkably similar to the pattern that the 
empirical correlation has tentatively singled out as the 
cause of compressor surge. This pattern is depicted in 
Fig. 14. In Fig. 15 the induced downwash from this vorticity 
pattern is pictorialized . If the blade tips were highly 
loaded, the additional lift caused by the circumferential 
vorticity definitely would be large enough to initiate stall. 

It seems, therefore, that there is a logical relationship be- 
tween a ring of large positive vorticity at the tip and stall 
that concurs with the purely empirical correlation which was 
pointed out in the previous section. For a thorough expla- 
nation of the method of calculation of these secondary effects, 
consult Ref. 8 or Ref. 1. 




FIG, M VORTICITY PATTERN 
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FIG. 15 INDUCED DOWNWASII [REF. 8] 




FIG. 16 REGIONS OF COMPRESSOR FACE 
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IV. A DISTORTION INDEX 



A. DEVELOPMENT OF THE INDEX 

From the two approaches to the interpretation, it has 
been seen that stall in the J-85 is probably caused by a 
type of distortion which is characterized by a large nega- 
tive vorticity at the midspan and a large negative vorticity 
at the tip. Both interpretations point to this conclusion, 
making it fairly certain that this path of investigation will 
bear fruit. In order to investigate the problem further and 
to determine if the conclusion that stall is caused by that 
one particular vorticity pattern is correct, it would be 
helpful to develop a distoartion index based on the assumption 
that that pattern causes stall. 

The distortion index will have certain advantages over 
the maps. The most important of these is that it will take 
less computer time and will facilitate the processing of more 
data. This will make it possible to evaluate data for which 
no map was produced and to check the conclusion made in the 
interpretation of the maps. 

The index should fulfill certain requirements. It should 
be large when the vorticity pattern, tentatively labeled as 
the stall producing pattern, is present and smaller when it is 
not. It should be large only if the pattern persists for one 
millisecond or more, since if it does not persist for that 
long, it is unlikely that it would have produced stall. 
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Since this is the first index ever devised based on vorticity, 
there is no guideline as to what is large and what is small. 
This is an advantage in that it does not matter how large the 
index is. It is only the relative size at different points 
in time that is of interest. However, it will be possible to 
get some feel for the size of this index which causes stall 
through the comparison of the index for different stall events 
The development of the index will now be described. 

Through inspection it was seen that the positive ring of 
vorticity was generally at r ' of greater than 0.75, and the 
negative ring of vorticity was at r 1 of less than 0.75. The 
boundary layer had dominated the vorticity at r' greater than 
0.95. At r 1 less than 0.5, there was very little regularity 
to the pattern. In order to reduce computer time, it was 
decided to employ the observations above and to make the in- 
dex a simple arithmetic summation with as little testing for 
magnitude as possible. To this end, the vorticity map was 
separated into three sections as shown in Fig. 16. Region I 
at r‘ less than 0.5 was ignored. In region II, the vorticity 
at radii of 0.55, 0.60, 0.65, and 0.70 was found for 0 of 0 
to 360 at steps of 0.5°, and summed yielding equation (4). 

7/7 

lAtnCt) i-Vcar^cKNU-s') r V^r (.cs^ks) rvk.’j0 7o ^ x-y j] (4) 

llzi 

If the vorticity at any point in this region was greater than 
0, it was rejected. In region III, the positive vorticity at 
radii of 0.8, 0.85, 0.90, and 0.95 was found and summed. In 
this case negative numbers were rejected. This yields 
equation ( 5 ) . 
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The distortion factor is then the difference between these 
two sums, (equation 6). 



At this point, a distortion factor was known at each 
instant of time. If the ring of positive vorticity at the 
tip is present along with the ring of negative vorticity at 
the hub, this index should be a large number. However, this 
number carries no information about the duration of the 
pattern. In this simple development the time factor was in- 
corporated in the following manner. The distortion factor 
for each case was multiplied by the distortion factor 0.625 ms. 
before it and 0.625 ms. after it. In this manner a new dis- 
tortion factor is calculated that reflects the vorticity pat- 
tern as seen through a "window" that is 1.25 milliseconds 
long. This is a simple matter to program since the cases used 
are all 0.625 milliseconds apart. The distortion factor 
becomes 



In this manner a particular distortion must last at least 
1.0 millisecond in order to yield a large value of the 
distortion factor. 




( 6 ) 
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B. TEST OF THE DISTORTION FACTOR 



The distortion factor was inserted into the computer 
program and tested for several stall events. A graphical 
representation of the results of each of the first four 
events is presented in Figs. 17-20. 

It can be seen that the distortion factor is quite 
effective for three of the four cases. In the fourth case, 
however, it gives no indication whatsoever of stall. This 
could be due to one or all of several influences. The index 
fails to account for several things. It could be that the 
index becomes large at the wrong time due to the fact that 
there may be a ring of positive vorticity at the tip and 
negative vorticity at the midspan without the increase in 
lift if they do not coincide as shown in Fig 21. If this 
is the case, they could cause the index to be large without 
inducing stall. Also, there may be other influences, that 
are at the present obscure, which play a part in the stall of 
a compressor. Finally, it could be that a mild distortion 
can cause stall if it persists for long enough. 

Another interesting point about the test is that for two 
of the three successful tests the peak of the distortion index 
which seems to have initiated stall is of the same order of 
magnitude, while the third is much smaller. This seeming 
discrepancy could provide some insight into the mechanism of 
stall. The operating points for the four cases were plotted 
on a performance map for the J-05, and the steady-state 



40 



operating point for the case in which a smaller peak in 
distortion index initiated stall is much closer to the surge 
line than the other two successful tests. This may be a 
significant observation. 

In any case, if the conclusions that lead to the dis- 
tortion index are correct, then the initial level of dis- 
tortion is not the same for all operating conditions of the 
compressor. Also, there is something about that particular 
operating point for the fourth case that makes the compressor 
more susceptible to stall. 




FIG, 21 VORTICITY PATTERN PRDUCES LARGE INDEX, 
NO STALL 
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FIG. 20 DISTORTION INDEX, CASE IV 



A 

I i 

I __ 
STALL— | 




CO 



V. SUMMARY AND REVIEW 



In review it is necessary to put the findings of this 
thesis in perspective. In order to test the usefulness of 
the instantaneous vorticity approach to the problem of in- 
let distortion in axial flow compressors, the pressure 
distribution at the compressor face must be converted to a 
vorticity distribution at the compressor face. A computer 
program was developed which performed this task with reason- 
able accuracy. 

Once the program was operational it was used to convert 
instantaneous pressure data for the J-85-GE-13, which was 
taken in the 10 x 10 foot supersonic wind tunnel at NASA 
Lewis Research Center in Cleveland, Ohio, into vorticity 
maps of radial and circumferential vorticity at the compres- 
sor face. The map of radial vorticity offered little infor- 
mation, but a marked empirical relationship was found between 
a circumferential vorticity pattern and stall. A crude 
instantaneous distortion index was defined using this empirical 
relationship as a criterion. The results of these tests show 
that there is indeed an empirical relationship between instan- 
taneous vorticity and stall. 

The tests of the distortion factor also showed that this 
pattern of vorticity is not the only mechanism that is capable 
of causing stall. 
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VI. CONCLUSIONS 



The primary purpose was to establish the validity of 
the vorticity approach to the problem of instantaneous dis- 
tortion in a jet inlet. This approach to the problem has 
proved to be quite useful in understanding the effects of 
distortion on the stability of a compressor. The computer 
program which converts the pressure distribution at the 
compressor face to the vorticity pattern at the compressor 
face will be instrumental in further research along this 
line of reasoning. 

While the crude distortion factor developed in this work 
is fairly successful in the prediction of stall, and could 
most likely be refined to be even more effective, it is not 
the purpose of this line of research to determine an empiri- 
cal instantaneous distortion pattern for the J-85. The data 
for the J-35 were used to show that the vorticity approach 
to the problem is valid. The overall purpose of the re- 
search is to develop a universal inlet distortion factor 
that could be employed for a wide variety of engines and 
airframes. Further research should be along these lines. 

It may be possible to make progress toward a universal 
inlet distortion factor through the calculation of the 
theoretical secondary flow effects of vorticity as it 
passes through the compressor. In further work the radial 
and axial vorticity must be incorporated into the analysis. 
Other factors that must be considered in the analysis include 
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blade twist and type of blading. Finally, as stated in 
Ref. 2, other types of distortion, such as temperature 
distortion, must be considered. The vorticity approach to 
the problem of inlet distortion is very promising and may 
provide the reliable universal distortion criteria that 
industry lacks today. 
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APPENDIX A 



EXPLANATION OF DATA 



The data used by the computer program are the same data 
which were used in Ref. 5. They were taken in the 10 x 10 
foot supersonic wind tunnel at NASA Lewis Research Center in 
Cleveland. They consist of eleven different stall events. 

For each stall event the 250 milliseconds of analog data 
prior to stall were digitized at 8000 points/second/channel. 

A schematic of the digitizing system is in Fig. A-l. Each 
of the digitized values consists of only the fluctuating 
component of pressure, so the steady-state values of pressure 
just prior to stall must be used in order to find the absolute 
total pressure at each point. The steady-state value was not 
used alone, however . Instead a reference pressure which in- 
cludes a correction for the leakage of the high pass filter 
was added to the fluctuating component to obtain absolute 
values of the stagnation pressure. The reference pressure 
was determined in the following manner. Analog data for all 
thirty channels were filtered at 50 hertz. A drawing of the 
compressor face with the positions of all the probes is shown 
in Fig. A-2. A five second' section of data prior to stall 
was digitized at 200 points/second/channel. The thousand 
values for each channel were then averaged to find delta p 
average for each channel, where delta p average is equal to 
the filter leakage. The values of P re f and dP/dt then were 
calculated using equations A-l and A-2. 
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FIG. A-l SCHEMATIC OF DIGITIZING SYSTEM 
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For each event the time of engine stall wa s known as well 
as the time of earliest indication of stall inside the 
compressor. This time of earliest indication of stall inside 
the compressor is the time at which the interstage dynamics 
indicated that the engine was going to stall and is accurate 
only to the nearest millisecond. A synopsis of these times 
for each stall event is in Table A-I. Also, a list of 
steady-state operating conditions for each stall event is 
contained in Table A-I I. 

If a more complete explanation of the instrumentation 
used is desired, consult Ref. 5. 
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APPENDIX B 



ACCURACY OF THE COMPUTER PROGRAM 



When a computer program of this complexity is used, it 
is mandatory that the accuracy of the program be verified. 

In the calculation of the vorticity distribution at the com- 
pressor face there are six quantities involved: r * , $ , 

normalized axial velocity, normalized total pressure, and the 
derivatives of the normalized pressure in the r* and theta 
directions. The distance r' is actually one of the defining 
values for any point on the compressor face, and gamma is 
merely the ratio of specific heats which is a well known func- 
tion of temperature. The four remaining quantities are cal- 
culated by the program as a function of r 1 and 9 that is pro- 
vided by the interpolation scheme. Since the normalized 
axial velocity is a function of two measured quantities, the 
average total pressure and the average static pressure, and 
the total pressure, it is obvious that if the interpolation 
scheme accurately computes the total pressure at each point, 
then the normalized axial velocity also will be correct. 

Thus we see that the accuracy of the entire program hinges 
on the accuracy of the interpolation scheme which calculates 
P‘ 5Lf , and 5j? . The interpolation scheme is a cubic spline 

> dv ' 1 d& 

fit in two variables. It fits a cubic polynomial to each 
three points in one dimension and then matches the slope of 
the curve at the points that are common to both polynomials. 
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There are two ways to test the accuracy of this inter- 
polation scheme. One way is to set up a dummy function of r' 
and Q, and produce a grid of point values of this dummy 
function. Then, using the interpolation scheme, predict the 
value of the dummy function and its two derivatives at other 
points in between the grid values and compare these values 
to the true values of the dummy function which are known 
since the dummy function is known explicitly. The second way 
is to plot the known values of the pressure on a graph, con- 
nect the points with a curve, and see if the interpolated 
values fall on the curve. The slope of the curve can be 
compared with the partial derivatives that were calculated. 
The results of the first method of testing are in Table B-l. 
In Figs. B-l - B-5 the results of the second method are 
presented . 

The interpolation scheme was found to have a maximum 
error of one percent for the value of pressure in this test, 
while the values of the partial derivatives were in error by 
a maximum of four percent. 

Once the accuracy of the interpolation scheme has been 
verified, the only remaining possibility of error is the sub- 
routines that print the map of the vorticity distribution. 

In Table B-II the values of vorticity calculated on a desk 
calculator are compared with the symbols that were printed 
at these points by the computer while drawing the maps. 

The positions of these points are indicated on the vorticity 
map in Fig. B-6. 
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TABLE B-I 



PREDICTED VALUES OF A FUNCTION VS ACTUAL VALUES 



Actual 


Predicted 


Actual 


Predicted 


Actual 


Predicted 


1918.2 


1920.5 


- 97.49 


-100.3 


45.66 


43.59 


1927.3 


1929.4 


-120.0 


-122.9 


141.5 


144 . 6 


1947.5 


1950.0 


-152.5 


-156.1 


268.6 


271.0 


1982.0 


1984 . 9 


-196.7 


-201.1 


426.8 


431.7 


2033.9 


2036.8 


-254.3 


-259.6 


616.2 


617.6 


2106.3 


2111.3 


-327.3 


-355.0 


836.8 


856.6 


1723.9 


1723.0 


-265.1 


-265.2 


-269.9 


-267.5 


1693.9 


1693.3 


-310.3 


-309.6 


-329.9 


-331.9 


1657.9 


1657.1 


-375.2 


-374 . 8 


-389.9 


-391.1 


1615.9 


1614 . 9 


-463.6 


-463.2 


-449.9 


-452.2 


1567.9 


1566. 8 


-578.9 


-578.1 


-509.9 


-511.6 


1513.9 


1512.1 


- 724.9 


-726.5 


-569.9 


-576.3 


1382.6 


1382.1 


-378.1 


-378.7 


-585.6 


-574.4 


1313.5 


1314.2 


-400.7 


-401.0 


-801.5 


-802.1 


1221.3 


1222.0 


-433.1 


-433.7 


-1048.6 


-1044.5 


1102.8 


1103.8 


-477.3 
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-1326.8 


-1324.4 
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-463.7 


-463.1 


-585.6 


-574.4 


872.7 


873.4 


-441.2 


-440.9 


-801.5 


-802.1 


780.5 


781.2 


-408.7 


-408.2 


-1048.6 


-1044 . 5 


661.9 
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-364.5 


-363.8 


-1326.8 


-1324 . 4 


514.1 


516.0 


-306.9 


-306.1 


-1636.2 


-1625.8 
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THETA 




FIG. B-l PRESSURE vs THETA, r'=OAI 
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FIG, B-2 PRESSURE vs THETA, r ' =0 . 59 
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FIG. B-3 PRESSURE vs THETA - r '=0 . 73 
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FIG. B-'l PRESSURE vs THETA- r '*=0 . 85 
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FIG, B-5 PRESSURE vs THETA, r'*0. 95 
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TABLE B-II 



TEST OF PRINT ROUTINE SYMBOLS 



r 1 0 W' Symbol Symbol Which 

Printed Should be Printed 
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